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Abstract — This paper considers a Popov type approach to the 
problem of robust stability for a class of uncertain linear quan- 
tum systems subject to unknown perturbations in the system 
Hamiltonian. A general stability result is given for a general 
class of perturbations to the system Hamiltonian. Then, the 
special case of a nominal linear quantum system is considered 
with quadratic perturbations to the system Hamiltonian. In this 
case, a robust stability condition is given in terms of a frequency 
domain condition which is of the same form as the standard 
Popov stability condition. 

I. Introduction 

This paper builds on the previous papers [l]-[3] which 
consider the problem of robust stability analysis for open 
quantum systems subject to perturbations in either the system 
Hamiltonian or coupling operator, which together define 
the dynamics of the quantum system. The results of these 
papers can be regarded as extensions of the classical small 
gain theorem for robust stability to the case of quantum 
systems. The main contribution of this paper is a result 
which can be regarded as an extension of the classical Popov 
criterion for absolute stability to the case of open quantum 
systems. In particular, we extend the result of [2], in which 
the perturbations to the system Hamiltonian are uncertain 
quadratic perturbations, to obtain a corresponding Popov 
robust stability result. 

The small gain theorem and the Popov criterion for 
absolute stability are two of the most useful tests for robust 
stability and nonlinear system stability; e.g., see [4]. Both 
of these stability tests consider a Lur'e system which is 
the feedback interconnection between a linear time invariant 
system and a sector bounded nonlinearity or uncertainty. 
The key distinction between the small gain theorem and the 
Popov criterion is that the small gain theorem establishes 
absolute stability via the use of a fixed quadratic Lyapunov 
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function whereas the Popov criterion relies on a Lyapunov 
function of the Lur'e Postnikov form which involves the sum 
of a quadratic term and a term dependent on the integral 
of the nonlinearity itself. The small gain theorem can be 
used to establish stability in the presence of time-varying 
uncertainties and nonlinearities whereas the Popov criterion 
only applies to static time-invariant nonlinearities. However, 
the Popov criterion is less conservative than the small gain 
theorem. Hence, we are motivated to obtain a quantum Popov 
stability criterion in order to obtain less conservative results. 

The study of quantum feedback control theory has been 
the subject of increasing interest in recent years; e.g., see [5]- 
[17]. In particular, the papers [14], [18] consider a framework 
of quantum systems defined in terms of a triple (S,L,H) 
where S is a scattering matrix, L is a vector of coupling 
operators and H is a Hamiltonian operator. The paper [18] 
then introduces notions of dissipativity and stability for this 
class of quantum systems. As in the papers [1], [3], the 
results of this paper build on the stability results of [18] to 
obtain robust stability results for uncertain quantum systems 
in which the quantum system Hamiltonian is decomposed as 
H = Hi + H-2 where H\ is a known nominal Hamiltonian 
and i?2 is a perturbation Hamiltonian, which is contained in 
a specified set of Hamiltonians W. 

For this general class of uncertain quantum systems, the 
paper first obtains a general abstract version of the Popov 
stability criterion which requires finding a Lyapunov type 
operator to satisfy an operator inequality. The paper then 
considers the case in which the nominal Hamiltonian H\ is 
a quadratic function of annihilation and creation operators 
and the coupling operator vector is a linear function of 
annihilation and creation operators. This case corresponds to 
a nominal linear quantum system; e.g., see [8], [9], [11], [12], 
[17]. Also, it is assumed that the perturbation Hamiltonian is 
quadratic but uncertain. In this special case, a robust stability 
stability criterion is obtained in terms of a frequency domain 
condition which takes the same form as the classical Popov 
stability criterion. 

The remainder of the paper proceeds as follows. In Section 
llll we define the general class of uncertain quantum systems 
under consideration. In this section, we also present a general 
Popov type stability result for this class of quantum systems. 
In Section [III] we consider a class of uncertain quadratic 
perturbation Hamiltonians. In Section IIVI we specialize to 
the case of linear nominal quantum systems and obtain a 
robust stability result for this case in which the stability 
condition is a frequency domain condition in the same form 
as the classical Popov stability condition. In Section [V] we 
present an illustrative example involving a quantum system 



arising from an optical parametric amplifier. In Section [VI] 
we present some conclusions. 

II. Quantum Systems 

We consider open quantum systems defined by parameters 
(S,L,H) where H = H l + H 2 ; e.g., see [14], [18]. The 
corresponding generator for this quantum system is given by 



G(X) = -i[X,H] + £(X) 



(1) 



where C(X) = \L*[X,L] + \[L*,X]L. Here, [X,H] = 
XH — HX denotes the commutator between two operators 
and the notation * denotes the adjoint of an operator. Also, 
H\ is a self-adjoint operator on the underlying Hilbert space 
referred to as the nominal Hamiltonian and H 2 is a self- 
adjoint operator on the underlying Hilbert space referred to 
as the perturbation Hamiltonian. The triple (S, L, H), along 
with the corresponding generators define the Heisenberg 
evolution X(t) of an operator X according to a quantum 
stochastic differential equation; e.g., see [18]. 

The problem under consideration involves establishing 
robust stability properties for an uncertain open quantum 
system for the case in which the perturbation Hamiltonian 
is contained in a given set VVi. Using the notation of [18], 
the set Wi defines a set of exosystems. This situation is 
illustrated in the block diagram shown in Figure Q] The main 



Exosystem 



Input Field 



Nominal 
System 



Output Field 



Fig. 1. Block diagram representation of an open quantum system interacting 
with an exosystem. 

robust stability results presented in this paper will build on 
the following result from [18]. 

Lemma 1 (See Lemma 3.4 of [18].): Consider an open 
quantum system defined by (S, L, H) and suppose there ex- 
ists a non-negative self-adjoint operator V on the underlying 
Hilbert space such that 



g(v) + cv < x 



(2) 



where c > and A are real numbers. Then for any plant 
state, we have 

(V(t)) <e- ct (V) + -, Vt >0. 

c 

Here V(t) denotes the Heisenberg evolution of the operator 
V and (■) denotes quantum expectation; e.g., see [18]. 



We will also use following result, which is a slight modifi- 
cation of Theorem 3.1 and Lemma 3.4 of [18]. 

Lemma 2: Consider an open quantum system defined by 
(S,L,H) and suppose there exists non-negative self-adjoint 
operators V and W on the underlying Hilbert space such 
that the following quantum dissipation inequality holds 

g(v) + w<\ (3) 

where A is a real number. Then for any plant state, we have 
1 



limsup 

T->oo 1 



(W(t))dt < A. 



(4) 



Here W(t) denotes the Heisenberg evolution of the operator 
W and (■) denotes quantum expectation; e.g., see [18]. 

Proof The proof is similar to the proof of Lemma 3.4 in 
[18]. In a similar manner we obtain from (|3) 



V(t + h)-V(t) + 



t+h 



W{s)ds 



< Xh. 



(5) 



where E t denotes vacuum expectation operator. Taking the 
vacuum expectation E on both sides, and noting that 
ErjE t = Eo results in 

(ip,E [V(t + h)]ip) - (ip,E [V(t)]tp) 

(•t+h 

(tp,E [W(s)]ip)ds < \h{i>, tjj) 

for any ip in the underlying Hilbert space. Then for any 
corresponding plant state (i.e., for ip such that = 1) 

we have 

±(V(t)) + (W(t))<\. 
Then follows in a standard manner. 



□ 



A. Commutator Decomposition 

We now consider a set of self-adjoint perturbation Hamil- 
tonians H 2 E Wi. For a given set of non-negative self- 
adjoint operators V, a set of Popov scaling parameters 
C [0, oo), a self-adjoint operator Hi, which is the nominal 
Hamiltonian, a coupling operator L, and for a real parameter 
P > 0, this set Wi is defined in terms of the commutator 
decompositions 

[V-9Hi,H 2 ] = [V -9Hi,7J]w- wt[z,V-6Hi], 
C{H 2 ) < £(z t )w + w t £(z) + f3[z,L}^[z,L] 

(6) 

for all V £ V and 9 £ 0, where w and z are given operator 
vectors of the same dimension. Here, the notation t denotes 
the adjoint transpose of a vector of operators. In addition, 
the notation # denotes the vector of adjoint operators for a 
given vector of operators. 

Then, the set Wi will be defined in terms of the sector 
bound condition 



1 

w ; 

7 



1 

w ; 

7 



< 



1 



z^ z 



(7) 



where 7 > is a given constant. That is, we define 

w f H 2 > : 3 w, z, such that \ 

1 \ and © are satisfied MV E V, 9 e 6 J ' ( ' 

Using this definition, we obtain the following theorem. 

Theorem 1: Consider a set of non-negative self-adjoint 
operators V, an open quantum system (S, L, H) and an 
observable W where H = H\ + H 2 and H 2 € Wi defined in 
dSJ. Suppose there exists a K G P and real constants 9 G 6, 
A > such that 

-»[V,fli]+£(V) 

+ i(i[z,l/-6»Hi]+6l/:(z) + z) t 

7 



Furthermore, 



x (i[ar, V - + &C(z) + z) 
+9/3[z,L] jf [z,L] + W < A. 



(9) 



Then 



limsup— / (W(t))dt < A. 

Here denotes the Heisenberg evolution of the operator 

W. 

Proof: Let 9 € 6 and A > be given such that the conditions 
of the theorem are satisfied and consider Q(V + 9H 2 ) defined 
in (Q~|i. Then 

g(V + 6H 2 ) = -i[V + 9H 2 ,Hi+H 2 }+£(V + 9H 2 ) 
= -i[V, Hi] - %B[H 2 ,Hx] - i[V, H 2 ] 

-t9[H 2 ,H 2 ]+C(V) + 9C(H 2 ) 
= -i[V,H!]-i[V -9H ll H 2 } + C(V) 
+9£{H 2 ). (10) 

Using the decomposition in the first equation (O, we have 

G(V + 9H 2 ) = -i^H^+L^-ilV -6H uZ t]w 
+iw*[z,V -9H 1 ] + 9C(H 2 ). (11) 

Now 

[V-QH^z^ = z(V-9H 1 )-(V-9H 1 )z = [z,V-BHi] 

since V — 9 Hi is self-adjoint. This confirms that the operator 
on the right hand side of the above identity is a self-adjoint 
operator. Therefore, the following inequality follows from 
the second equation ©: 

Q{V + 9H 2 ) 
< -i\V,Hi] + C(V) 

-i[V - 9H llZ ^]w + iw 1 [z, V - 6Hx] 
+6>(£(z t )w + w*£(z)) + 9/3[z, L}*[z, L]. (12) 

Also, note that 

(£(z))t = (\tf[z,L\ + \[L\z\L) 



(13) 



= l[L\zl]L + \tf[z\L\ 
= £(zt). 







x 

1 

1 



i[z,V -9H 1 ]+9C(z) 
i[z,V -9H l ]+9C{z 



(i[z,V -9H 1 ]+9C(z)+z) 



y/l(w- zli) 
t 



x {i[z, V -6H X ] +9C(z) + z) 
- (i[z, V - 6H X ] + 9C{z) + z) f w 
-w f (i[z, V - 9H X ] + 9C{z) + z)+ jw^w. 

This implies that 

V - 0.ffi] +9C(z) + z) t w 
W (i[z, V - 9H{\ + 9C(z) +z)- jw^w 

< -(i[z,V-9H 1 ]+9C(z) + z) 1 

7 

x {i[z,V ~9H 1 ]+9C(z) + z). 
Using this inequality and (fT2l i. we have 

G(V + 9H 2 ) + --y(w-^z 

< -i[V,Hi)+£(V) 

+ (i[z, V - OHx] + 9C(z) + z) f w 
W (i[z, V - 9H X \ + 9C(z) + z) 
—jw^w + 6/3[z,L]^[z,L] 

< -i\V.H 1 ]+ C{V) 

+ -(i[z,V ~9H 1 \+9C{z) + z) ] 

7 

x (i[z,V -6H 1 ) + 6C(z) + z) 
+9j3[z,L]^[z,L}. (14) 

Then it follows from (|7]l and (O that 

g(V + 9H 2 ) + w<\. 

The result of the theorem then follows from Lemma [2] □ 

III. Quadratic Perturbations of the Hamiltonian 

We consider a set W 2 of quadratic perturbation Hamilto- 
nians of the form 



Ho 



A 



c 



(15) 



where A € 



^2mX2m 



A 



is a Hermitian matrix of the form 

" Ai A 2 " 
Af Af 



(16) 



Ai = Aj, A2 = A 2 and £ is a given vector of operators. 
Here, in the case of complex matrices, the notation t refers 
to the complex conjugate transpose of a matrix. Also, in 
the case of complex matrices, the notation # refers to the 
complex conjugate matrix. In addition, for this case we 
assume that 9 = [0, 00). 



The matrix A is subject to the bounds 

4 

< A < -I. 

7 

Then we define 

H 2 of the form dTSb . (IT6b such that 



W 2 



condition ( fTTI i is satisfied 



(17) 



(18) 



Using this definition, we obtain the following lemma. 

Lemma 3: Suppose that [z, L] is a constant vector. Then, 
for any set of self-adjoint operators V, 



W 2 c Wi. 

Proof: Given any H 2 E W 2 , let z — 

Hence, H 2 = w^z. Then, for any V E V and 8 > 0, let 
V = V — 9Hi and we have 



and w = \Az. 



[V,z*}w = -(Vz^-z^V^Az 

= -Vz^Az - -z^AVz, 
2 2 

since V is a scalar operator and A is a constant matrix. Also, 



■[z,V] = i*t A (. 



V- Vz 



-z^AzV — -z^AVz. 



Hence, 



[V,z f )w- w f [z,V] 



Ivz^Az-lz^AzV 



= \[V,^Az] 
= [V,H 2 \. 

Similarly, 

[L,z*]w-w*[z,L] = [L,H 2 ]. 

In addition, 

w [z,L]*L-L[z,L]t w 
= \z T A T [z,L]#L-\L[z,L^A; 

= -[z,L]^AzL- -L[z,L]^ALz 



(19) 



and similarly 

L*[z,Lfw* 



,L]^A[z,L] 
,L]*A[z,L] 



,L}L* 



(20) 



1 



z,L]*A[z,L] 



Now using j!9\ . d20l ). (fSTJ and the assumption that [z, L] is 
a constant vector, it follows that 

C{H 2 ) = ^L*[H 2 ,L] + ^[L*,H 2 }L 



L* ([L, 



z'\w — w 



--( W i[L,z^-[z,L?w)L 
= ±L'[zi,L]w + ±L'w*[z,L] 

+-wt[z\ltf L + -[z,ltfwL 
= \l*[z\L]w + \l*[z,L] t w* 

+ ly[z\L}^L + ^w T [z,L]*L 

= ^L*[z\L]w + ^[z,L}L* 

+±w1[z\L]iL + ±L[z,Lpw 
+ [z,L] t A[z,L] 
= h*[z\L]w + ^[L*,zJ]Lw 

+ ^L*[z,L] + ^[L*,z}L 

+[z,L?A[z,L] 
= C(zi)w + w*£(z) + [z,L]tA[z,L]. 
It then follows from (TTtT > that 

C{H 2 ) < C(z^)w + w^C(z) + -[z,Lp[z,L]. 

7 

Therefore we can conclude that both of the conditions in (|6]i 
are satisfied with (3 = —. Also, condition ( fTTI i implies 



H 2 = wh = -z^Az > 0, 



and 



1 y / 1 

w z w z 

7 / V 7 



t 1 t 
w 1 w z 1 w 

7 

1 t 1 t 
w' z H — kZ' z 

7 7 

-z^AAz - -z^Az 

4 7 

1 t 

H — ^z'-z 
7 



< 



-z' z 



r 



which implies (0. Hence, F 2 G W\. Therefore, W2 C Wi. 

□ 

IV. The Linear Case 

We now consider the case in which the nominal quantum 
system corresponds to a linear quantum system; e.g., see [8], 
[9], [11], [12], [17]. In this case, we assume that H\ is of 
the form 



(21) 



Fx = 1 [ ffl t a T ] M 



a 



(22) 



where M £ 



••2nx2n 



is a Hermitian matrix of the form 



M 



M 1 M- 



AI* M- 



1-2 
T* 



and Mi = M[, M 2 = Mj. Here a is a vector of annihilation 
operators on the underlying Hilbert space and a# is the 
corresponding vector of creation operators. The annihilation 
and creation operators are assumed to satisfy the canonical 
commutation relations: 























a 




a 




a 




a 




. a * . 




. a * 






. a * . 




. a * . 





a 


# 


a 




( 


. a * . 




. a * . 


1 



= J 



(23) 



where J 







-I 

In addition, we assume L is of the form 



e.g., see [10], [15], [17]. 



L = [ Ni N 2 
where Ni G C lxn and N 2 G 



JV 



(24) 





= N 


a 






a* 





Also, we write 



Ni N 2 



N* NI 



it 





a 




a* 



In addition, we assume that V is of the form 

V = [ at a T ] P 



a 



(25) 



where P G 
of the form 



-i2rix2n 



is a positive-definite Hermitian matrix 



P = 



Pi 



Pi 

5 # 



(26) 



Pi 



off 
2 ^1 

Hence, we consider the set of non-negative self-adjoint 
operators Pi defined as 

V of the form (f25]l such that P > is a 
Hermitian matrix of the form (l26l 

'(27) 

In the linear case, we also let £ = E\a + E 2 a# and hence 
we can write 

Ei E 2 
Ef Ef _ 

(28) 

We will also consider a specific notion of robust mean 
square stability. 

Definition 1: An uncertain open quantum system defined 
by (S, L, H) where H = Hi + H 2 with Hi of the form 
(l22l . H 2 G W, and L of the form (l24l is said to be robustly 
mean square stable if for any H 2 G W, there exist constants 
Ci > 0, c 2 > and C3 > such that 



c 











a 


= E 


a 




. a * . 







a(t) 


t 


a(t) 


( 


a#{t) 




a#(t) 



< cie 



a 



a 



c 3 Vt > 0. (29) 



Here 



a(f) 
a#(t) 



denotes the Heisenberg evolution of the 



vector of operators 



a 

n* 



; e.g., see [18]. 



In order to address the issue of robust mean square 
stability for the uncertain linear quantum systems under 
consideration, we first require some algebraic identities. 

Lemma 4: Given V E V\, Hi defined as in (l22l i and L 
defined as in d24l i. then 



[V,Hi 



a 

a* 



a 

n* 



Also, 

C(V) 



[PJM - MJP] 
1 



1 

' 2 



M 



a 

n* 



a 

a* 



-IS[V,L] + -[L\V]L 



Tr PJN^ 



I 




N J 

(NUNJP + PJN^JN) 



Proof: The proof of these identities follows via straightfor- 
ward but tedious calculations using d23l . □ 
Lemma 5: With the variable z defined as in d28l i and L 
defined as in d24"l i. then 



*,L] 



E 



a 

n* 



,N 



a 



= EJUN 



which is a constant vector. Here, 



E = 



I 

1 



Proof: The proof of this result follows via straightforward 
but tedious calculations using ( f23l . 

Lemma 6: With the variable z defined as in 
defined in ( f22b and L defined as in d24l , then 



□ 
Hi 



-i[z,Hi]+C(z) 



where 



Furthermore, 



= E 
= EA 



-UM 
a 



2 J 



a 

a* 



UM - -JN^JN. 



(30) 



i[z,V]=2iEJP a# 

Proof: The proofs of these equations follows via straightfor- 
ward but tedious calculations using d23l . □ 
We will show that a sufficient condition for robust mean 
square stability when H 2 G W 2 is the existence of a constant 
9 > 0, such that the following conditions are satisfied: 

(i) The matrix A defined in ( f30b is Hurwitz. 



(ii) The transfer function 

G(s) = —2iE (si - A)" 1 JE^ (31) 
satisfies the strict positive real (SPR) condition 

7/ - (1 + 6iu)G(iu>) - (1 - 6iu)G{iurf > (32) 

for all cj G [—00, 00]. 

This leads to the following theorem. 

Theorem 2: Consider an uncertain open quantum system 
defined by (S,L,H) such that H = Hi + H 2 where Hi is 
of the form (|22]i, L is of the form d34j and H 2 G W 2 . 
Furthermore, assume that there exist a constant > 
such that the matrix A defined in (f30b is Hurwitz and 
the frequency domain condition (l32l is satisfied. Then the 
uncertain quantum system is robustly mean square stable. 
Proof of Theorem [2] If the conditions of the theorem are 
satisfied, then the transfer function ZI — (1 + 6s)G(s) is 
strictly positive real. However, this transfer function has a 
state space realization 



1 



jI-(l + 6s)G(s) 



where A is defined as in d30l ), 

B = -2iJE^ 

and 

C = E. 



A 


B 


-c - ecA 


\i - ecB 



(33) 



It now follows using the strict positive real lemma that the 
linear matrix inequality 



< 



PA + A^P PB + C^ + 6A^C^ 

B^P + C + OCA -7/ + 6{CB + B t C t ) 

(35) 

will have a solution P > of the form ( f26b ; e.g., see [4]. 
This matrix P defines a corresponding operator V € Vi as in 
( |25"1 >. Furthermore, it is straightforward to verify that CB + 
B^C^ = 0. Hence, using Schur complements, it follows from 
that 



PA + A^P 

+-(PB 

7 



Now using Lemma [6] we have 

i[z,V -9Hi] + 6£(z) + z 
= i[z,V]+6{-i[z,Hi] + C(z)) + z 



(2iEJP + 9EA + E) 
E {2UP + 6A + I) 



a 

a* 



a 

a* 



Hence using Lemma [H we obtain 



-i[V,Hi]+£(V) 
1 



6 Hi 



+ -(i[z,V 

7 

x (i[z,V-6Hi\ 
AG 

+-[z,L]1[z,L] 
7 



+ ec{z) + z) ] 
ec(z) + z) 



a 

n* 



M 



Tr f PJN f 
46» 



,# 



/ 





N*T 1 JE^EJY,N T 



NJ 



1 



(37) 



where 

M = 

PA + A^P + 

- (2iJP + 6A + J) t E^E (2iJP + 6A + I) 

7 

= PA + A^P 

+- (PB + C 1 + OAW) (B*P + C + 9CA) , 

A is defined in (f30b . B is defined in d3~3l and C is defined in 
d34l i. From this, it follows using Lemma [5] Lemma [3] d36l l. 
and a similar argument to the proof of Theorem Q] that 

t 



g(v + 0H 2 ) < 



a 

a* 



M 



(I 

a* 



X 



(34) where 



Tr ( P JN 1 



I 




NJ 



+—N*Y<JE^EJY.N T 

7 

> 0. 

It follows from d36t that M < 0. Hence using ( fTTb , it follows 
that there exists a constant c > such that the condition 

G(V + 6H 2 ) + c{V + 6H 2 ) < X 

is satisfied. Therefore, it follows from Lemma Q] Lemma [3] 
<H2) and P > that 



» ( 


o(t) 


f r (t) 




)« 






I a # w J 


(36) 




o(0) 


t 


" o(0) 






a#(0) _ 




_ a#(0) 



10 



E^E] 



^min [P] 



cX 

nil it 

[P] 



w > 0. 



(38) 



Hence, the condition $2% is satisfied with ci 

x min \p\ 1 > 0, c 2 = c> and c 3 = cA J i Pj > 0. □ 

Observation 1: A useful special case of the above result 
occurs when the QSDEs describing the nominal open quan- 
tum linear system depend only on annihilation operators and 



not on the creation operators; e.g., see [11], [12]. This case and 
corresponds to the case of M2 = and N2 = 0. Also, we 
assume that E2 = 0. In this case, we calculate the matrix A 
in (O to be 

" A 1 
A* 

where A x = -«Mi-|iV|iVi. Also, we calculate the transfer Also A 



.4 



function matrix G(s) in OTb to be 



G(s) 



Gi(s) 






-Gi(a* 



where G x (s) = iE x {si - A x )~ x El. 

We now consider the case in which G x {s) is a SISO 
transfer function. In this case, the condition that the matrix 
A in d30b is Hurwitz reduces to the condition that the matrix 

1 



A x 



-iMi 



(39) 



is Hurwitz. Also, the the SPR condition d32l reduces to the 
following conditions: 



Z+He[Gi(iu)] 



6u)Tm[Gi(iui)] 
6(dTm[Gi(iu>)] 



> 



> 



(40) 
(41) 



for all uj G [—00,00]. The conditions d40l . (|4TT > can be 
tested graphically producing a plot of uj7Ze[Gi(iuj)] versus 
lm[Gi(iuj)] with uj G [—00,00] as a parameter. Such a 
parametric plot is referred to as the Popov plot; e.g., see 
[4]. Then, the conditions d40l . d41~h will be satisfied if and 
only if the Popov plot lies between two straight lines of slope 
i and with x-axis intercepts see Figure [2] 



allowable region 



uiTm[G\{iuj)\ 




slope = I 

7ee[Gi(iw)] 



Fig. 2. Allowable region for the Popov plot. 

V. Illustrative Example 

In this section, we consider an example of an open 
quantum system with 

S = I, H = —i ^(a*) 2 — a 2 ^ , L = y/na, 

which corresponds an optical parametric amplifier; see [19]. 
This is the same example which was considered in [2]. In 
order apply the theory of this paper to this example, we let 



Ho = ha" 





1 i ' 




a 


a] 


-i 1 




a* 



> 



Hi = -[a* a] 



" -1 




a 


-1 




a* 



so that H1 + H2 = H. This defines a linear quantum system 
of the form considered in Theorem |2] with M\ = —1, 

M 2 = 0, JVi = y/n, N 2 = 0, Ex = 1, E 2 = 0. 
1 i 

> 0, which satisfies A < 21. 



-i 1 j 

Hence, we can choose 7 = 2 to ensure that condition (TTtT > 
is satisfied and therefore H2 G W^- Also, note that this 
system is a system of the form considered in Observation 
Q] with A x = i — #, which is Hurwitz for all k > 0, and 
Gi(s) = We then choose k = 2.1 and construct the 

S 2+ 2 

Popov plot corresponding to the transfer function Gi(s) as 
discussed in Observation Q] For a value of = 0.2, this plot, 
along with the corresponding allowable region, is shown in 
Figure [3] From this figure it can be seen that the Popov plot 
lies within the allowable region and hence, it follows from 
Theorem [2] and Observation Q] that this system will be mean 
square stable for k = 2.1. However, the method of [2] could 
only prove that this system was stable for k > 4. Hence for 
this example, the proposed method provides a considerable 
improvement over the method of [2]. 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 



Fig. 3. Popov plot for the optical parametric amplifier system. 



VI. Conclusions 

In this paper, we have considered the problem of robust 
stability for uncertain linear quantum systems with uncertain 
quadratic perturbations to the system Hamiltonian. The sta- 
bility condition which is obtained is a quantum version of 
the classical Popov stability criterion. This frequency domain 
condition is less conservative than a previous stability result 
obtained which takes the form of a quantum version of the 
classical small gain theorem. 
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